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Abstract 

Using 281 pb“^ of data collected with the CLEO-c detector, we report on first observations and 
new measurements of Cabibbo-suppressed decays of D mesons to 2, 3, 4, and 5 pions. Branching 
fractions of previously unobserved modes are measured to be: B{D^ 7r''“7r“7r^7r°) = (9.9 ± 
0.6 ± 0.7 ± 0.2 ± 0.1) X 10-3, tt+tt+tt-tt-tt^) = (4.1 ± 0.5 ± 0.2 ± 0.1 ± 0.0) x 10-3, 

B{D+ tt+ttOttO) = (4.8 ± 0.3 ± 0.3 ± 0.2) x 10-3, B{D+ 11 + 71 + 71 - 71 ^) = (11.6 ± 0.4 ± 0.6 ± 
0.4) X 10-3, B{D^ giT^) = (0.62 ± 0.14 ± 0.05 ± 0.01 ± 0.01) x 10-3, and B{D^ ujtt+tt-) = 
(1.7± 0.5 ± 0.2 ± 0.0 ± 0.0) X 10-3. The uncertainties are from statistics, experimental systematics, 
normalization and CP correlations (for modes only). Improvements in other multi-pion decay 
modes are also presented. The D tttt rates allow us to extract the ratio of isospin amplitudes 
A{AI = 3/2)/A{AI = 1/2) = 0.420 ± 0.014(stat) ± 0.016(syst) and the strong phase shift of 
Si = (86.4 ± 2.8 ± 3.3)°, which is quite large and now more precisely determined. 
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Charm decays provide an important laboratory for the stndy of both the strong and weak 
interactions. Semileptonic decays provide direct access to CKM matrix elements. Hadronic 
decays provide important inpnt to B physics as well opening a window into the stndy of 
hnal state (strong) interactions. Exhanstive or precise measnrements of Cabibbo-snppressed 
decays have been challenging dne to low rates or other experimental challenges, however, 
more information on these decays is of great importance in several areas. Firstly, in the 
weak sector, Cabibbo-snppressed hnal states, snch as —>• 7 r+ 7 r“ 7 r°, provide a promising 
way by which to extract the CKM angle 7 in B'^ —> ^pg similar magnitnde 

of interfering amplitndes between and > Tr+Tr^vr® In the arena of 

strong interactions, they provide a means in which to stndy hnal state interactions [ 3 , , 

which may help clarify how they contribnte to i? —> tttt decay |^. Moreover these modes can 
enter as backgronnds to other D decay measnrements and therefore their valnes are of 
general importance to ascertain. Information on these decays is sparse and not very precise. 
The CLEO-c '0(3770) ^ DD sample provides the opportnnity to perform a comprehensive 
and precise stndy of these decays. 

The CLEO-c detector is a general pnrpose solenoidal detector which inclndes a tracking 
system for measnring momenta and specihc ionization of charged particles, a Ring Imaging 
Cherenkov detector to aid in particle identihcation, and a Csl calorimeter for detection of 
electromagnetic showers. The CLEO-c detector is described in detail elsewhere j^. 

This analysis ntilizes 281 pb“^ of data collected on the 0(3770) resonance (-^/s ~ 3.77 
GeV) at the Cornell Electron Storage Ring. At this energy, DD pairs are prodnced in a 
coherent 1 hnal state with no additional particles. 

We reconstrnct 77° mesons in several mnlti-pion decay channels, inclnding D'^ , 

D^ —> Tr'^TT^, D^ —TT+TT^TT^, D^ —Tr+TT+TT^Tr”, D^ —> Tr+TT^TT^Tr^, n~'K ~and 

D^ —> In D^ we reconstrnct D^ —>■ D^ D^ —>■ D^ —>■ 

7 r+ 7 r’'' 7 r“ 7 r°, D^ —>■ 7 r+ 7 r+ 7 r+ 7 r“ 7 r“. We also consider the resonant p and uj contribntions and 
measnre rates for D^ —> D^ —»• a; 7 r''' 7 r“, and D^ r]7T~^. The branching fractions are 

measnred relative to D^ —> K~7i~^ and D~^ —>■ for nentral and charged D mesons, 

respectively. Thronghont the paper, charge conjngate modes are implicitly assnmed, nnless 
otherwise noted. 

The reconstrnction of D mesons uses charged particles and vr^’s reconstructed with stan¬ 
dard selection requirements [3|. For each candidate, we utilize two kinematic variables: 

AE = Ebeam “ Ed and M^c = \JEl^^^ - Pd, where Ed (pd) are the energy (momentum) 
of the D candidate and Ebeam is the beam energy. The substitution of the beam energy for 
the candidate energy improves the mass resolution by about a factor of hve. For properly 
reconstructed candidates, AE exhibits a narrow peak near zero and Mbc peaks at the D 
meson mass. Mode-dependent signal regions in AE are dehned to be within about three 
times the r.m.s. widths of the AE distribution obtained from detailed Monte Carlo (MC) 
simulations (see Tabled). A AE sideband region extending from 10 MeV beyond the signal 
region to 100 MeV is used to study the shape of the background in Mbc. 

The Cabibbo-favored modes such as D —y Kg -|- X, Kg —>■ tttt constitute a large back¬ 
ground as their branching fractions are typically hve to ten times larger than the Cabibbo- 
snppressed modes, and therefore such events must be vetoed. For each decay channel with 
> three pions, we veto any candidate which contains a pair of pions with invariant mass in 
the range 475 < < 520 MeV/c^ or 448 < M^o^o < 548 MeV/c^. This range was de¬ 

termined from a large sample of generic MC events where we require that the surviving Kg 
background be less than 1% of the expected signal in the corresponding Cabibbo-snppressed 
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FIG. 1: Mbc distributions for modes from data. The points are the data and the superimposed 
lines are the fits as described in the text. 


signal channel. 

The resulting Mbc distributions for the neutral and charged D modes are shown in Figs.^ 
and|2l respectively. The points show the data from the /S.E signal region and the lines are hts 
to the distributions which are given by the sum of an ARGUS threshold function and an 
asymmetric signal shape(CBAL) which models the initial state radiation (ISR) induced 
tail at large Mbc- The ARGUS shape parameters are extracted by htting the /S.E sidebands. 
The signal shape parameters are determined by hts to Mbc distributions obtained from a 
simulation of DD production at the ■^(3770) jlO| followed by a simulation of the detector 
response mill- The Gaussian widths of the Mbc distributions, which are part of the 
GBAL signal shape, and the htted yields are listed in Table Q] The yields we hnd for the 
normalization modes are in good agreement with the previously published values 0. 

Efficiencies for a given decay mode depend mildly on the presence of intermediate reso¬ 
nances, but the dependence is increased to as large as 10% (relative) when the veto is 
included. To minimize this bias, we tune the MG simulation by either using existing Dalitz- 
plot analyses, as for the TrTrvr hnal states umi, or by adding intermediate resonances to 
each mode in order to reproduce the observed substructure in the tttt invariant mass distri¬ 
bution. In most cases, this requires introducing signihcant p contributions. Efficiencies from 
simulation are checked in data by comparing the numbers of particles of each species in fully 
reconstructed events with the corresponding yields when their populations are inferred from 
energy and momentum conservation [^. The resulting corrections are: (-3.9±2.0)% per 7r°, 
(-0.7±0.7)% per and (-0.3±0.3)% per [3|. For vr^’s with momentum near Md/2, 
the correction is (—4.6 ± 3.5)%. The resulting efficiencies are shown in the last column of 
Table m 

These multi-pion hnal states are fed by intermediate resonances, such as p, rj, cn, 0, /o, 
02, etc. We examine the multipion hnal states for p, a; —only. The 4 


TT 


TT 71 ^ 
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FIG. 2: Mbc distributions for D~^ modes from data. The points are the data and the superimposed 
lines are the fits as described in the text. 


TABLE I: Requirements on AE for signal candidates, Gaussian widths of the Mbc distributions, 
observed yields, and reconstruction efficiencies. See text for details. 


Mode 


AE 

Tfbc 


Efficiency 



[low, high] Width 

Yield 

(%) 



(MeV) 

(MeV) 



DO ^ 

TT~^TT 

[-30, 30] 

1.42 

2085±54 

73.1±0.7 

DO ^ 

O 

o 

[-50, 40] 

2.88 

499±32 

31.0±0.7 

DO ^ 

+ 

1 

o 

[-50, 40] 

1.83 

10834±164 

39.9±0.7 

DO ^ 

+ 

+ 

1 

1 

[-25, 25] 

1.46 

7331±130 

48.4±0.6 

DO ^ 

+ 

1 

o 

o 

[-40, 30] 

2.02 

2724±166 

15.9±0.5 

DO ^ 

+ 

+ 

1 

1 

o 

[-30, 20] 

1.67 

1614±171 

18.4±0.6 

DO ^ 

o 

o 

o 

[-60, 35] 

3.00 

29±15 

13.4±0.7 

DO ^ 

K~7r~^ 

[-29, 29] 

1.42 

51210±231 

65.0±0.6 

D+ ^ 

> 7r'*“7r0 

[-40, 35] 

1.98 

914±46 

46.5±0.8 

D+ ^ 

> TT^TT^Tr” 

[-25, 25] 

1.40 

3303±95 

62.0±0.8 

D+ ^ 

+ 

o 

o 

[-30, 30] 

1.99 

1535±89 

22.0±0.6 

D+ ^ 

• TT+TT+TT TfO 

[-30, 30] 

1.58 

5701±205 

30.6±0.7 

D+ ^ 

+ 

+ 

+ 

1 

1 

[-15, 15] 

1.38 

732±77 

27.6±0.6 

D+ ^ 

> K 7r'’“7r“’“ 

[-22, 22] 

1.37 

80381±290 

54.4±0.5 


intermediate state is not treated here since it can be measured significantly better using 
0 —>■ K~^K~, and the observed rates are too low to improve on existing measurements. 

We search the —>■ 7T~^7T~^n~7r^ , —>■ Tr+Tr^Tr^TT®, and 7T~^'n'^n~7r~7T^ modes for rj 

and uj decays. The yields are extracted by selecting events that are within 2.5 times the r.m.s. 
width of the D mass and taking the difference in yields between the number of such events 
in the AE signal and AE sideband regions. The sideband distributions are normalized to 
account for the different range of AE between signal and sidebands regions. Distributions of 
7 T~^7r~7T^ mass for these three modes are shown in Fig.|21 The left column shows combinations 
whose AE is in the signal region and the right column shows candidates in the AE sideband 
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FIG. 3: Distributions in 7r'*'7r“7r® invariant mass for —> 7r'’“7r'’“7r“7r’^ (a-jb), —> 7r'’'7r“7r^7r’^ 

(c,d), and —> 7r'*“7r'’'7r“7r“7r*^ (e,f) candidates. The left column represents the Aif signal region 

and the right are for I\E sidebands. The superimposed lines are the fits as described in the text. 


region. The distributions are fit to the sum of a polynomial background and two Gaussians, 
one each for the rj and u. With the limited statistics and large background, we do not try 
to £t the Breit-Wigner tails of the ut, but rather assume a Gaussian shape and absorb the 
loss of events from the tails into the efficiency. In the £t, the rj and uj Gaussian widths are 
constrained to 3.8 MeV and 9.0 MeV, respectively, the values determined from signal MG. 
The resulting yields for the AE signal and sideband regions are shown in Table O along 
with the reconstruction efficiency, determined from signal MG. 


TABLE II: Summary of submodes containing r] and u) mesons, indicating efficiencies and yields in 
the AE signal and AE sideband regions. 


Mode 


Efficiency 

Yield 

Yield 



(%) 

AE signal AE sideband 

D+ 

► r]7r~^ 

29.4±1.0 

421±23 

44T12 

D+ ^ 

► 

21.7±1.0 

216±43 

236±41 


rjTT^ 

21.9±1.0 

90±12 

28±8 


UJTT^ 

14.1±1.0 

103±26 

140±25 


1 

+ 

20.3±1.0 

260±32 

150±29 


(UTT+Vr 

15.6±1.1 

1304±96 

832±91 


Relative branching fractions are computed using the yields shown in Tables Q] and El 
and are presented in Table |TT ] funseen decay modes of the r] and oj are included, using 
B{ri,u —>• n~^n~Tr^) from Ref. 15[). To compute the absolute branching fractions, we use 
B{D^ —>■ K~7r~^) = (3.84±0.07)% and ^ iT“7r+7r+) = (9.4±0.3)%, which are obtained 
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using a weighted average of the PDG values and the recent CLEO measurements 0]. For 
unobserved modes, we set 90% conhdence level upper limits. In the last column of Table HTTl 
we show PDG ^ averages, when available. 


TABLE III: Measured relative and absolute branching fractions for neutral and charged D modes. 
Uncertainties are statistical, experimental systematic, normalization mode uncertainty, and un¬ 
certainty from CP correlations (for modes only). For the relative branching fractions, the 
normalization mode uncertainty is omitted. 


Mode 


^mode/^ref 

^mode 

B (PDG) 



(%) 

(10-3) 

(10-3) 

7)0 ^ 

TT+vr 

3.62 ±0.10 ±0.07 ±0.04 

1.39 ± 0.04 ± 0.04 ± 0.03 ± 0.01 

1.38±0.05 

7)0 ^ 

O 

o 

2.05 ±0.13 ±0.16 ±0.02 

0.79 ± 0.05 ± 0.06 ± 0.01 ± 0.01 

0.84±0.22 

7)0 ^ 

+ 

1 

o 

34.4 ±0.5 ±1.2 ±0.3 

13.2 ±0.2 ±0.5 ±0.2 ±0.1 

11±4 

7)0 ^ 

+ 

+ 

1 

1 

19.1 ±0.4 ±0.6 ±0.2 

7.3 ±0.1 ±0.3 ±0.1 ±0.1 

7.3±0.5 

7)0 ^ 

+ 

1 

o 

o 

25.8 ±1.5 ±1.8 ±0.3 

9.9 ±0.6 ±0.7 ±0.2 ±0.1 


7)0 ^ 

o 

1 

1 

+ 

+ 

10.7 ±1.2 ±0.5 ±0.1 

4.1 ±0.5 ±0.2 ±0.1 ±0.0 


7)0 ^ 

LOTT^TT 

4.1 ± 1.2 ±0.4 ±0.0 

1.7 ±0.5 ±0.2 ±0.0 ±0.0 


7)0 ^ 

rjTT^ 

1.47 ±0.34 ±0.11 ±0.01 

0.62 ± 0.14 ± 0.05 ± 0.01 ± 0.01 


7)0 ^ 

O 

O 

o 

- 

< 0.35 (90% CL) 


7)0 ^ 

wttO 

- 

< 0.26 (90% CL) 


7)0 ^ 

1 

+ 

- 

< 1.9 (90% CL) 


7)+ - 

> TT^VrO 

1.33±0.07±0.06 

1.25±0.06±0.07±0.04 

1.33±0.22 

7)+ - 

> 1 T^TT^TT~ 

3.52±0.11±0.12 

3.35±0.10±0.16±0.12 

3.1±0.4 

7)+ - 

+ 

o 

o 

5.0±0.3±0.3 

4.8±0.3±0.3±0.2 


7)+ - 

> TT^TT^TT TfO 

12.4±0.5±0.6 

11.6±0.4±0.6±0.4 


7)+ ^ 

+ 

+ 

+ 

1 

1 

1.73±0.20±0.17 

1.60±0.18±0.16±0.06 

1.73±0.23 

7)+ ^ 

> r/7r+ 

3.81±0.26±0.21 

3.61±0.25±0.23±0.12 

3.0±0.6 

7)+ - 

> iOTr~^ 

- 

< 0.34 (90% CL) 



Several sources of systematic uncertainty are considered. Where appropriate, we include 
in parentheses the minimum to maximum size of the uncertainty. Limited MG statistics are 
used in determining the reconstruction efficiencies, which introduce relative uncertainties at 
the level of (1-3)%. Tracking and vr® reconstruction efficiencies have relative uncertainties of 
0.7% and 2.0% (3.5% for high momentum 7r°) per particle, respectively, and charged particle 
identihcation is understood at the level 0.3%/7r^ and 1.3%/iP^. Uncertainty from the Kg 
veto is estimated as the difference in probabilities (between data and MG) for each hnal state 
to pass the Kg veto. This survival fraction, is given by = (1 — /veto)^'"“L where 
/veto represents the veto probability per tttt pair, which is obtained by a linear interpolation 
from an 80 MeV region above and below the veto region, and iVpair is the number of tt+tt” (or 
7r°7r°) pairs in the given decay mode. The uncertainties on the veto efficiencies for 7r’''7r“ and 
7r°7r° are added in quadrature when applicable (0.0%-2.8%). Uncertainties in signal yields 
are estimated by comparing changes in the branching fraction when (a) signal widths are 
permitted to float, (b) varying the ISR-tail shape parameters individually by ±1 standard 
deviation, and (c) varying the range of the ht to the Mbc distribution. The three sources are 
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added in quadrature to obtain the total uncertainties (1.0%-4.9%). Uncertainties resulting 
from different AE resolutions are quantified by widening the AE windows, recomputing the 
branching fractions, and taking the difference between the new and default values ( 0 . 2 %- 
7.7%). The simulation of hnal state radiation has been studied in —>• events and 

is estimated to introduce no more than 0.5% uncertainty on the D reconstruction efficiency. 
Possible bias due to resonant substructure is estimated by removing the Kg veto (so as not 
to double count this systematic) and comparing the efficiencies obtained using our default 
simulation (which includes resonances) and a phase space simulation (0.0%-2.2%). Imperfect 
replication of the average number of candidates per event by the simulation could lead to 
a bias in the reconstruction efficiency. This effect is quantihed by comparing the average 
number of D candidates per event, (A^cand), within 2.5 standard deviations of the D mass, 
between data and simulation (on a mode-by-mode basis). The systematic uncertainty is 
taken as the difference in the ratio from unity, i.e., — 1 (0.2%-2.7%). The 

and normalization modes introduce uncertainties of 1.8% and 3.2%, respectively. Lastly, 
the effects of quantum correlations of DD pairs produced from the 1 "0(3770) may shift 

the branching fractions with respect to the values in the absence of these correlations. The 
effect of these quantum correlations has been considered Hi. and a shift in the branching 
fraction can be expected iiy = AT/T % 0. Limits on y are at the percent level 0 , and we 
take this as an additional systematic uncertainty on the branching fractions. The total 
systematic uncertainties on the branching fractions for each mode are shown in Table IIIIl 

Using the new H tttt branching fractions and D lifetimes, T£,+ = (1040±7) fs and = 
(410.3 ± 1.5) fs 113 , we compute the ratio of the AI = 3/2 to AI = 1/2 isospin amplitudes 
and their relative strong phase difference 0 to be A 2 /AQ = 0.420 ± 0.014(stat) ± O.Ol(syst) 
and cos 5 1 = 0.062 ± 0.048(stat) ±0.058(syst). The large phase shift, <5/ = (86.4 ±2.8 ±3.3)°, 
indicates that hnal state interactions are important in H —> tttt transitions. These results 
represent a considerable improvement over previous measurements from CLEO Q and 
FOCUS 0 , both of which are consistent with our data. 

In summary, we report new measurements of six Cabibbo-suppressed decay branching 
fractions of D mesons and eight additional measurements, of which all except for > 7 r’'' 7 r“ 
and —>■ 7 r’'' 7 r’'' 7 r’'' 7 r“ 7 r“ provide large improvements over the existing world average values. 
The large value of the strong phase shift, 5/ obtained in the D —> tttt decay supports the 
conclusion that hnal state interactions are important in D decays. 
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the National Science Foundation, and the U.S. Department of Energy. 
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